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Abstract 

We experimentally investigate the depolarizing power and the polarization entropy of a broad 
class of scattering optical media. By means of polarization tomography, these quantities are 
derived from an effective Mueller matrix, which is introduced through a formal description 
of the multi-mode detection scheme we use, as recently proposed by Aiello and Woerdman 
( arXiv:quant-ph/0407234| ). This proposal emphasized an intriguing universality in the polarization 
aspects of classical as well as quantum light scattering; in this contribution we demonstrate exper- 
imentally that this universality is obeyed by a surprisingly wide class of depolarizing media. This, 
in turn, provides the experimentalist with a useful characterization of the polarization properties 
of any scattering media, as well as a universal criterion for the validity of the measured data. 

PACS numbers: 42.25.Dd, 42.25.Ja, 42.81.-i 



I. INTRODUCTION 



Characterization of optically transparent media with polarized-light scattering methods, 
is important for communication technology, industrial and medical applications [l^. When 
polarized light is incident on an optically random medium it suffers multiple scattering and, 
as a result, it may emerge partly or completely depolarized. The amount of depolarization 
can be quantified by calculating either the entropy {Ep) or the degree of polarization (Pp) 
of the scattered field jj]. It is simple to show that the field quantities Ep and Pp are related 
by a single- valued function: Ep{Pp). For example polarized light (Pp = 1) has Ep = 
while partially polarized light (0 < Pp < 1) has 1 > Ep > 0. When the incident beam is 
polarized and the output beam is partially polarized, the medium is said to be depolarizing. 
An average measure of the depolarizing power of the medium is given by the so called index 
of depolarization (-Dm) Non-depolarizing media are characterized hj Dm = 1, while 
depolarizing media have < Z^a/ < 1. A depolarizing scattering process is always accompa- 
nied by an increase of the entropy of the light, the increase being due to the interaction of 
the field with the medium. An average measure of the entropy that a given random medium 
can add to the entropy of the incident light beam, is given by the polarization entropy Em 
^. Non-depolarizing media are characterized by Em = 0, while for depolarizing media 
< Em < 1. As the field quantities Ep and Pp are related to each other, so are the medium 
quantities Em and Dm- In a previous paper 5] we showed the existence of a universal 
relation Em (Dm) between the polarization entropy Em and the index of depolarization Dm 
valid for any scattering medium. More specifically, Em is related to Dm by a multi-valued 
function, which covers the complete regime from zero to total depolarization. This universal 
relation provides a simple characterization of the polarization properties of any medium, as 
well as a consistency check for the experimentally measured Mueller matrices. We emphasize 
that the results found in apply both to classical and quantum scattering processes, and 
might therefore become relevant for quantum communication optical applications, where 
depolarization is due to the loss of quantum coherence 0. 

In this contribution, we present an experimental study of the depolarizing properties of a 
large set of scattering media, ranging from milk to multi-mode optical fibers. The results 
confirm the theoretical predictions for the bounds of the multi-valued function Ej^j{Dm)- 
The manuscript is divided as follows: in Section II we review the Mueller-Stokes formalism 
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and show the differences between deterministic (non-depolarizing) and non-deterministic 
(depolarizing) scattering media. We also discuss the statistical nature of a depolarizing 
process resulting from the average (either spatial or temporal) performed in a multi-mode 
detection scheme. Furthermore, in order to describe the transverse spatial average present in 
our multi-mode detection set-up, we formally introduce the concept of an effective Mueller 
matrix (Mes). In Section III we describe the experimental scheme for polarization tomog- 
raphy that was used to characterize the different scattering samples. These can be divided 
into two categories: (a) non-stationary (samples which fluctuate during the measurement 
time) and (b) stationary (samples which do not fluctuate). We then show the experimental 
results obtained for these samples followed by a brief discussion of the interesting structures 
in the {Em, Dj^j) plane, that were revealed by the experiments. Finally, in Section IV we 
draw our conclusions. 

II. DEPOLARIZING AND NON-DEPOLARIZING MEDIA 

In the Introduction we stressed the fact that passive optical systems may be grouped 
in two broad classes: depolarizing and non-depolarizing systems. To the first class belong 
all media which decrease the degree of polarization Pp of the impinging light, while to the 
second one belong all media which do not decrease Pp- In this Section we want to make the 
discussion more quantitative by using the Mueller-Stokes formalism which is widely used for 
the description of the polarization state of light beams. 

A. Mueller-Stokes formalism 

Consider a quasi-monochromatic beam of light of mean angular frequency lj Q]. Let 
us denote with x, y, z the axes of a Cartesian coordinate system, with the ^-axis along the 
direction of propagation of the beam whose angular spread around z is assumed to be small 
enough to satisfy the paraxial approximation. Let 



be the component of the complex paraxial electric field vector in the x- and y-direction 
respectively, at the point (x, y) located in the transverse plane z = zq at time to- If the field 




(1) 
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is uniform on the transverse plane, then Ex and Ey will be, in fact, independent of x and y 
and a complete description of the field can be achieved in terms of a doublets E of complex 
variables (with possibly stochastic fluctuations): 




where Eq and Ei are now complex-valued functions of Zq and to only. A complete study of 



the propagation of E along z can be found, e.g., in [8|, however, for our purposes the main 
result we need is that propagation through non- depolarizing media can be described by a 
deterministic Mueller (or Mueller- Jones) matrix M"^ , while to describe the propagation of a 
light beam through a depolarizing medium it is necessary to use a non- deterministic Mueller 
matrix M. 

1. Deterministic Mueller matrix M'^ 

In a wide-sense, a deterministic linear scatterer as, e.g., a quarter- wave plate, a rotator 
or a polarizer, is an optical system which can be described by a 2 x 2 complex Jones P] 
matrix 

^00 Joi 




(3) 

With this we mean that if E and E' describe the polarization state of the field immediately 
before and immediately after the scatterer respectively, then they are linearly related by the 
matrix J: 

E' = JE. (4) 

An alternative description can be given in terms of the Stokes parameters of the beam. To 
this end let C be the covariance matrix of the field defined as 

Q, = {E,E*), (z,j=0,l), (5) 

where the brackets denote the statistical average over different realizations of the random 
fluctuations of the field. Then the four Stokes parameters 5*^ (/x = 0, . . . , 3) of the beam are 
defined as 

S, = TT{Ca^}, (/i = 0,...,3), (6) 
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were the symbol Tr{-} denote the trace operation and the {cT/j} are the normahzed Pauh 
matrices: 





(7) 

i r)' 

Now, if with 5"^ and 5"^ we denotes the Stokes parameters of the beam before and after the 
scatterer respectively, it is easy to show that that they are linearly related by the real- valued 
4x4 Mueller- Jones matrix M"^ as 



where summation on repeated indices is understood and 



(8) 



(9) 



where the symbol "(8)" denotes the outer matrix product and the unitary matrix A is defined 
as 
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(10) 



From the structure of M"^ follows that a deterministic medium does not depolarize, that is 
Pf{S) — Pf{S') where the degree of polarization Pp of the field is defined as 

^Sf + SI + SI 



Pf{S) 



So 



(11) 



Let us conclude by noticing that for deterministic media the two descriptions in terms of J 
or M'^ arc completely equivalent in the sense that the 16 real elements of M'^ do not contain 
more information than the 4 complex elements of J. 



2. N on- deterministic Mueller matrix M 



A non-deterministic scatterer is, in a wide-sense, an optical systems which cannot be 
described by a Mueller- Jones matrix. In this class fall all the depolarizing optical system as, 
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e.g., multi-mode optical fibers, particles suspensions, etc.. It has been shown 11 1 that 
it is possible to describe a non-deterministic optical system as an ensemble of deterministic 
systems, in such a way that each realization £ in the ensemble is characterized by a well- 
defined Jones matrix J{£) occurring with probability p£ > 0. Then, the Mueller matrix M 
of the system can be written as 



(12) 



where the bar symbol denotes the average with respect to the ensemble representing the 
medium: 



At this point it is useful to introduce the auxiliary 4x4 Hermitian matrix H defined as 



H^„y = JijJli, (/i = 2i + j, iy = 2k + l), 



(14) 



which is, by definition, positive semidefinite, that is all its eigenvalues {Aq, Ai, A2, A3} are non- 
negative. Then, it is possible to show that the depolarization index Dm and the polarization 
entropy Em can be written as 



D 



M 



u=0 



1/2 



(15) 



^A/ = - J]A,log4(A,). (16) 

n 

From ref. |5| we know that Em is a multi- valued function of Dm and that this depen- 
dence determines some physical bounds to any polarization scattering process. The func- 
tion Em{Dm) shows thus character of universality. In the next Section we shall confirm this 
theoretical prediction with experimental results. 



B. Unpolarized light and depolarizing media 

In classical optics, a light beam appears to be depolarized when its polarization direction 
varies rapidly with respect to other degrees of freedom that are not measured during the 
experiment (e.g. wavelength, time or position of the beam) [l^. Moreover, depolarization 
occurs also when a single-mode polarization input beam is coupled with a multi-mode (either 
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spectral or spatial) system as, e.g., an optical fiber. In fact it is possible to identify two basic 
depolarizing processes, (a) one intrinsic to the medium, and (b) one due to the measurement 
scheme. In the first case (a) non- stationary temporal fluctuations of the optical properties 
of the medium, for instance due to the Brownian motion of suspended particles in a liquid 
[isj ]. cause depolarization even when a single- mode detection scheme is employed (the time 
average performed during the measurement is responsible for the depolarization). On the 
other hand, type (b) stationary depolarizers (i.e. glass fibers) do not fluctuate in time and 
produce light depolarization only in the presence of a multi-mode detection scheme. In this 
case it is simple to build explicitly the ensemble of Mueller- Jones matrices representative of 
the medium, that we introduced in the previous subsection. 

To this end, let us consider the case of a scattering process in which a coupling between 
polarization and spatial modes of the field occurs and a multi-mode detection scheme is em- 
ployed. This is, in fact, the case occurring for the optical fibers we used. The Mueller- Stokes 
formalism, is suitable for a single-mode description of the field; however, it is possible to 
extend this formalism to the case in which Ni^ spatial modes of the field impinge on the 
scatterer, A^^out leave from it and D modes are eventually detected. We make the assump- 
tion that different spatial modes of the field are uncorrelated, that is we do not consider 
interference phenomena which are ruled out by the required linearity with respect to the 
intensities of the field Moreover, without loss of generality, we assume N^^ = Nont = ^■ 
Let S(j) = {So{j), 5*1 (j), S2{j), 5*3 (j)} be a generic 4-D Stokes vector defined with respect to 
the mode j, where j G {1, . . . , A^}. For a A^-mode field we have a collection of of these 4-D 
Stokes vectors that we can arrange in a single 4A^-D "super" vector S = {S(1),..., S(A^)}. 
When the A^-mode light beam undergoes a polarization-sensitive scattering, then, in general, 
the Stokes vectors {Sin(j)} of the input beam are related to the set of vectors {Sout(j)} of 
the output beam by: 

N 

SoutU) = Yl Jo)S.„(io), (j = 1, • • • , N), (17) 

io=i 

where M'^(j,jo) is the 4x4 Mueller- Jones matrix that describes the scattering from the 
input mode jo to the output mode j. If we introduce a "super" 4A^ x 4A^ Mueller matrix 
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defined as 



( M-^(l,l) ... M\\,N)^ 



M = 



(18) 



\m-\N,V) ... M\N,N)^ 
where each block M"^(ji, J2) is a 4 x 4 Mueller- Jones matrix, then we can rewrite Eq. ipTjl 
in a compact form as 

^ ■ S^„. (19) 



'out 



After the scattering process took place, we have to detect its products. We recently showed 
that when D < N modes of the field are detected, a mode-insensitive polarization 
analyzer, put in front of a bucket-detector, can be described by a AN x AN block-diagonal 
matrix A: 



A = 



A{1) 



\ 



A{D) 



\ 







(20) 



where A{j), {j = 1, . . . ,D) are 4x4 real- valued positive semi-definite matrices (in fact, 
projectors), and is a null (A^ — D) x [N — D) matrix. In the paraxial limit {D << N) 
each A{j) reduces to the 4x4 identity. Then, the polarization state of the scattered beam 
after the analyzer, is described by the super Stokes vector S/j given by 



§D = A • M ■ Sin. 



(21) 



Finally, because of the mode-insensitive detection, the sum over all the detected modes 
reduces the number of degrees of freedom of the field from 4A^ to 4, producing the detected 
4-D Stokes vector 

N 

SD = Y,^DU)=M,ffS,n{l), (22) 

i=i 

where we have assumed that the input light beam is prepared in the single mode jo = 1, so 
that Sj„(jo) = SjgiSin{l) and with Meg we have denoted an effective 4x4 Mueller matrix 
defined as 

D 

" ' (23) 



which is written as a sum of D Mueller- Jones matrices. It is important to notice that while 
the product of Mueller- Jones matrices is still a Mueller- Jones matrix (in physical terms: 
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PU AU 

FIG. 1: Schematic of the polarization tomography set-up. For details see text. 

a cascade of non-depolarizing optical elements is still a non-depolarizing optical system), a 
sum, in general, is not. This causes depolarization. Moreover, since the "matrix coefficients" 
A{j) are non-negative, the matrix M^s in Eq. ()23|) is an explicit version of the Mueller 
matrix written in Eq. (fT^ . Then we have shown, by an explicit derivation, how to build 
the statistical ensemble representing the depolarizing medium, for this particular case. 

III. DEPOLARIZATION EXPERIMENTS 

A. Experimental scheme for polarization tomography 

In order to measure the effective Mueller matrix M^s and thus the index of depolarization 
Dm and the entropy Em of a scattering medium, it is straightforward to follow a tomography 
procedure: The light to be scattered by the sample is successively prepared in the four 
polarization basis states of linear {V, H, +45°) and circular (RHC) polarization, which are 
represented by four independent input Stokes vectors Sin. For each of these input fields the 
corresponding Stokes vector Sout, that represents the output field, is obtained by measuring 
the intensities of the scattered light in the same four polarization basis states. This procedure 
provides the 4x4 independent parameters required to determine the 16 elements M^j, of 
the Mueller matrix from Eq.(jH)). Note that we actually employ two additional polarization 
basis states (—45°, LHC) in our experiments and perform 6x6 measurement, which allows 
us to reduce experimental errors by averaging within the over-complete data set. 

The experimental scheme is illustrated in Fig. ^ The light source is a power-stabilized 
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He-Ne laser at 633 nm wavelength. The input field is prepared by the polarizer unit (PU), 
consisting of a fixed polarizer (PI), a half- wave plate (HI), and a quarter- wave plate (Ql). 
A microscope objective (MO, x 50/0.55) couples the fight into the sample. The scattered 
light is collimated by a standard photographic objective (PO, 50 mm/1.9), followed by an 
adjustable pinhole (PH) that defines the amount of transverse spatial average to be per- 
formed in the light detection. The analyzer unit (AU) consists of a quarter-wave plate (Q2) 
and a polarizer (P2). Together with a focusing lens (L) and a photodiode (PD), it probes 
the polarization state of the scattered output field. As an estimation of the systematic error 
of the set-up, mainly due to imperfections of the used retarders, we measured the Mueller 
matrix of air (i.e. the identity matrix) and of well-known deterministic optical elements such 
as wave-plates. In all these cases, we found the deviations from the theoretically predicted 
matrix elements limited to |AM^,^| < ±0.04. 

B. Collection of scattering media 

The various scattering media we investigated can be divided into (a) non- stationary 
samples where, e.g., Brownian motion induces temporal fiuctuations within the detection 

integration time, and (b) stationary samples without such fluctuations, most notably multi- 
mode polymer and glass optical fibers. More specifically, we chose our scatterers from: 

(a) Non-stationary media: 

— polystyrene microspheres (2 /im dia., suspended in water, Duke Scientific Co., 
USA). 

— diluted milk; 

(b) Stationary media: 

— Zenith™ polymer sheet diffusers (100 /im thick, SphereOptics Hoffman GmbH, 
Germany) ; 

— holographic fight shaping diffusers (0.5°, 1°, 5°, and 10° scattering angle. Physical 
Optics Co., USA); 

— quartz/silica wedge depolarizers and quartz Lyot depolarizers, Halho Optics, UK); 
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— step-index polymer optical fiber (NA=0.55, core dias. 250 /im, 500 /im, 750 /xm 
ESKA CK type, Mitsubishi Rayon, Japan); 

— step-index glass optical fiber (NA=0.48, core dias. 200 /im, 400 /xm, 600 /xm 
FT-x-URT type, distributed by Thorlahs, Inc., USA); 

— step-index glass optical fiber (NA=0.22, core dia. 50 /xm, ASF50 type, distributed 

by Thorlahs, Inc., USA); 

— graded-index glass optical fiber (NA=0.27, core dia. 62, 5 /im, GIF625 type, 
distributed by Thorlahs, Inc., USA). 

C. Experimental results 

For a large collection of different samples. Fig. |21 shows the measured polarization entropy 
Em vs. the corresponding index of depolarization Dm- The black lines represent the calcu- 
lated analytical boundaries in the {Em, Dm) plane, whose functional dependence Em{Dm) 
was derived in Ref. [Sj. These boundaries provide universal constraints to the possible values 
{Em, Dm) for any physical scattering system, that is, the range of admissible values is re- 
stricted to the rather limited grey-shaded area within the boundaries. As it is apparent from 
the experimental data, our choice of samples allowed us to widely fill in the range of values 
{Em, Dm), in good agreement with the prediction from Ref. For rather different scatter- 
ing media, we observed similar values of the pairs {Em, Dm), which display the universality 
in this quantitative description of the depolarizing properties. We found samples throughout 
the full range of values in entropy and depolarizing power, < Em, Dm < 1- However, note 
that the region below the curve connecting the points A and C in the {Em, Dm) plane is 
not covered by any data so far. Work is in progress to investigate this peculiarity. 
Several scatterer-specific tuning parameters allowed us to realize this wide range of depolar- 
izing systems and to reveal details of the depolarizing properties for the various media, as 
will be discussed in the following subsection. 

The most versatile scatterers used to acquire data in the {Em, Dm) plane, were the multi- 
mode optical fibers. For them, the depolarization is caused by multiple reflections within 
the cylindrical light-guide together with mode mixing. We selected fibers of various lengths 
that displayed the full range between non- depolarizing and completely depolarizing prop- 
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FIG. 2: Measured entropy (Em) vs. index of depolarization (Dm) for (a) non-stationary and 
(b) stationary scattering media. The maximal possible parameter range in the {Dm, Em) plane is 
indicated by the grey-shaded area. Lines correspond to analytical boundaries predicted by theory. 
Cuspidal points are given by A = (0, 1), B = (l/3,log4 3), C = (l/\/3, 1/2), and D = (1,0) fl]. 

erties. Fibers shorter than about 2 cm showed negligible depolarization {Dm ~ 1). In 
the case of the glass fibers complete depolarization {Dm ~ 0) was observed for lengths of 
5 m, whereas in the case of the polymer fibers this was achieved already for lengths of 
only ^ 50 cm. The reason is, presumably, the significant Rayleigh scattering at density 
fiuctuations in the polymer material [isj ]. 

In our experimental scheme, the aperture of the pinhole (PH) defines the region of spatial 
of average in the scattered light detection (see Fig. . By choosing the pinhole diameter 
between 2 mm and 13 mm, we realized scattering systems which are described by different 



12 




FIG. 3: Measured entropy {Em) vs. index of depolarization (Dm) for (a) step-index polymer 
optical fibers (♦); (b) step-index glass optical fibers (•), graded-index glass optical fibers (o); 
(c) microspheres (a), milk (t), polymer diffuser (□), holographic diffusers with wave plate (■), 
wedge depolarizers (a), and Lyot depolarizers (v). The analytical boundaries are indicated by the 
continuous lines. 

effective Mueller matrices Mes- In fact, a small pinhole, corresponding to an average over a 
small set of modes j in Eq. (jSHl), leads to a large index of depolarization. However, due to 
the huge optical mode volume in our fibers, the pinhole adjustment allowed only for small 
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changes, within |AZ}m| < ±0.05. With a lower hmit of 2 mm in the pinhole diameter, 
special care was taken to select a sufficiently large number of speckles in the scattered out- 
put field. (The step-index glass fiber with 50 /im core showed the largest speckles of about 
500 /im FWHM.) This is necessary in order to average out interference effects, generated 
by the coherent source (He-Ne laser), so that the assumption of uncorrelated modes in the 
derivation of M^s holds 3]. 



D. Discussion 



In Fig. 01 we separately show the results for (a) step-index polymer fibers, (b) step-index 
and graded-index glass fibers, and (c) other scattering media. For the polymer fibers it 
is apparent that most of the data fall on the upper curve connecting the points A and 
D. This curve corresponds to Mueller matrices that have an associated operator H with 
its four eigenvalues of the degenerate form {A, /x, /i, /i}. This degeneracy can be associated 
with isotropic depolarizers, which is obviously a good description for the polymer fibers. 
Contrarily, in Fig. Efb) (glass fibers), we ffil in also the allowed {Em, Dm) domain below 
the isotropy curve. These domains correspond to anisotropic media, the anisotropy being 
supposedly due to stress-induced birefringence in the glass fibers. It was actually this bire- 
fringence which we used as an additional tuning parameter accounting for changes of a few 
percent in the index of depolarization. The data obtained for long fiber samples, both in 
polymer and glass, are close to the cuspidal point A {\ = fi = 1/4), which corresponds to 
total depolarizers. Contrarily, the data for very short fiber samples are close to the cuspidal 
point D (A = 1, /i = 0), which corresponds to a deterministic (i.e. non-depolarizing) system. 
In case of watery suspensions such as milk and microspheres, we observed purely isotropic 
depolarization, i.e. all data are found on the isotropic curve, similar to the polymer optical 
fibers (see Fig. Efc)). We adjusted the depolarizing power by varying the concentration of 
the scatterers. The 100 fim thick Zenith"""^ polymer diffuser sheet was found to be almost 
completely depolarizing, whereas the holographic diffusers, when used alone, did not cause 
any significant depolarization. The latter effect is due to the absence of multiple scattering 
in the transmission of light through these surface-optical elements. In combination with a 
subsequent wave-plate we could, however, couple the scattered light angular spectrum to 
the polarization degrees of freedom, and thus achieve depolarization. 
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Finally, the data for standard wedge and Lyot depolarizers are also shown in Fig. Efc). 
Wedge depolarizers are designed to completely depolarize a well defined linear input po- 
larization, whereas our tomographic measurement procedure represents an average over all 
independent input polarizations. This results in a non-zero index of depolarization. The 
Lyot depolarizer is, within the experimental error, non-depolarizing since it is designed to 
depolarize a broadband light source while we operated with a monochromatic laser source. 



IV. CONCLUSIONS 



By means of polarization tomography we have experimentally characterized the depolar- 
izing properties for a large set of scattering optical media. We describe these media with 
both the index of depolarization Dj^i and the polarization entropy Em that is added by 
the medium to the scattered field. These quantities are derived from a measured effective 
Mueller matrix, which we formally introduce in the description of scattering systems 
subject to a multi-mode detection, as is the case of our experimental configuration. The 
set of studied media ranges from non-stationary scatterers such as milk and polystyrene 
microspheres to stationary scatterers such as multi-mode optical fibers, diffusers, and 
standard wedge depolarizers. 

In Ref. 0] a universality was predicted for the possible values of Em and Dm, these values 
being restricted to a limited domain described by a set of analytical boundaries. The 
collected experimental data for our scatterers fill in this domain almost completely and give 
evidence of the predicted depolarization universality in light scattering. A certain range of 
the predicted {Em, Dm) values is, however, not covered by the scatterers we investigated so 
far. Work on this issue is currently under progress. 

Furthermore, the quantities Dm and Em provide insights into the particular depolarization 
mechanisms of the various media, as well as a consistency check for the measured data (see 
), and may provide a useful classification of optical scatterers for quantum applications, 
where depolarization stands for decoherence j?]. In this spirit, an extension to twin photon 
quantum scattering experiments is relatively straight-forward. Work along this line is also 
under progress in our group. 
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